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The momentum distribution of quark-gluon plasma at the early stage of a relativistic heavy- 
ion collision is anisotropic and consequently the system, which is assumed to be weakly coupled, 
is unstable due to chromomagnetic plasma modes. We consider a high-energy parton which flies 
across such an unstable plasma, and the energy transfer between the parton and the medium is 
studied as an initial value problem. In the case of equilibrium plasmas, the well-known formula 
of collisional energy loss is reproduced. The unstable plasma case is much more complex, and the 
parton can lose or gain energy depending on the initial conditions. The extremely prolate and 
extremely oblate systems are considered as examples of unstable plasmas, and two classes of initial 
conditions are discussed. When the initial chromodynamic field is uncorrelated with the color state 
of the parton, it typically looses energy, and the magnitude of the energy loss is comparable to that 
in an equilibrium plasma of the same density. When the initial chromodynamic field is induced by 
the parton, it can be either accelerated or decelerated depending on the relative phase factor. With 
a correlated initial condition, the energy transfer grows exponentially in time and its magnitude 
can much exceed the absolute value of energy loss in an equilibrium plasma. The energy transfer 
is also strongly directionally dependent. Consequences of our findings for the phenomenology of jet 
quenching in relativistic heavy-ion collisions are briefly discussed. 


I. INTRODUCTION 


The observation of jet quenching in central collisions of relativistic heavy ions is considered to be a signal that 
quark-gluon plasma (QGP) is produced, because only a medium with deconfined color charges could stop a multi- 
GeV parton within a few femtometers, see e.g. the reviews mm- The energy loss of a high energy parton is a key 
element of a quantitative understanding of the jet quenching phenomenon and has been intensively studied for two 
decades, but the problem is far from being completely solved mm- 

QGP is produced in the early stage of a relativistic heavy-ion collision and spends most of its lifetime in a state 
of local equilibrium. The energy loss of a high-energy parton is therefore conventionally computed in a locally 
equilibrated plasma which evolves hydrodynamically. However, a QGP reaches a state of local equilibrium only after 
a short but finite time interval mm, and the momentum distribution of pre-equilibrium plasma is anisotropic. Even 
if the anistropic phase is very short lived, it might have a significant effect on the energy loss of a test parton because 
the weakly coupled anisotropic QGP is unstable due to spontaneously growing chromomagnetic modes (for a review 
see [5]). The anisotropic plasma is populated with large chromodynamic fields which will strongly influence the test 
parton. Therefore, the energy loss of the parton in the brief anisotropic phase might constitute a significant fraction 
of the total energy loss responsible for the experimentally observed jet quenching. An analysis of this system is very 
complicated however, because the unstable QGP evolves quickly due to the presence of unstable modes, and the 
energy loss calculation has to be treated as an initial value problem. Our results show that the energy loss is strongly 
time dependent, and this dependence is much stronger than the switching-on effect studied in mm- In contrast to 
the weakly coupled plasma studied here, energy loss in the strong coupling regime computed in the framework of 
AdS/GET duality is rather similar in equilibrium and in far-from-equilibrium plasma at the same energy density [H] . 

In this paper we study a highly energetic parton which gains or loses energy through interactions with the chro- 
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modynamic fields present in the QGP. The effect of elastic interactions is called collisional energy loss and that of 
gluon emission radiative energy loss. In the case of light partons - light quarks and gluons - radiative energy loss is 
expected to give the dominant contribution to the total energy loss mm- For heavy quarks the radiative energy loss 
is presumably less important due to the so-called dead cone effect mm- Collisional energy loss in anisotropic QGP 
was studied in Ref. [3] but the unstable plasma was treated as a static medium and the interaction of the test parton 
with exponentially growing chromodynamic fields was ignored. The energy loss computed in [9] thus misses the key 
feature of unstable plasmas. 

We will treat the test parton as a classical particle with SU(A^c) color charge, the dynamics of which is described by 
the Wong equations m combined with the linearized Yang-Mills equations. We assume that the typical momenta of 
the collective modes are much less than the typical momenta of the plasma constituents. This approach is equivalent to 
using QCD within the hard-loop (HL) approximation [11] . In the equilibrium limit the time dependence of the energy 
loss disappears and we reproduce the soft part of the collisional energy loss [I2HIS], where the momentum transfer 
is of the order of the Debye mass. The energy loss due to soft interactions diverges logarithmically with the upper 
limit of the momentum transfer, which we call /cmax- In our anisotropic calculations, we also find an approximately 
logarithmic dependence on fcmax- This ultraviolet sensitivity is expected since the approach is classical, and signals the 
necessity to combine the classical contribution to the energy loss at small wave vectors with the quantum contribution 
at higher ones. A quantum approach to parton energy loss in unstable plasma needs to be developed. 

Our crucial finding is that depending on the initial conditions the test parton can either lose or gain energy when 
traversing the unstable QGP. In an equilibrium plasma the parton loses energy and the energy change per unit length 
dE/dx is negative. If the parton gains energy from the plasma fields, dE/dx is positive. Although the energy transfer 
can be either negative or positive, depending on the situation, we frequently use the term ‘energy loss’ generically to 
describe both cases. Our results show that when the intial conditions are chosen in a certain way, the magnitude of the 
energy loss increases exponentially, which indicates that the unstable modes play an important role. At late enough 
times, the energy loss can be much bigger than in equilibrium plasmas. It is also strongly directionally dependent. 

The acceleration of a test particle in a plasma system might seem rather exotic, but the phenomenon is well known 
in the physics of electromagnetic plasmas. It is caused by the electric field associated with plasma waves in the system. 
Gharged particles are carried forward on the electrostatic wave with a motion like surfing with speed equal to the 
phase velocity of the wave, and can therefore be boosted to very high energies. This picture motivates the idea to use 
a plasma excited by a laser or particle beam as a particle accelerator. A mechanism was proposed in 1979 by Tajima 
and Dawson m. and was experimentally verified soon afterwards m Since plasmas can sustain accelerating fields 
orders of magnitude larger than those in the radio-frequency modules of standard accelerators, small plasma devices 
can be extremely efficient. In the experiment in Ref. [18], electrons were accelerated to an energy as high as 1 GeV 
over a distance of 3.3 cm, demonstrating immense promise for affordable and compact plasma accelerators for various 
applications. 


This paper is organized as follows. In Sec. |n| we derive the general energy loss formula of a relativistic classical 
parton in an unstable QCD plasma which depends on the initial conditions. The method is similar to the approach 
developed earlier by one of us to study the spectra of chromodynamic fluctuations [19] and the momentum broadening 
of a fast parton [^. The equilibrium limit is discussed in Sec. Ill where we show that the dependence on the initial 
conditions drops out, and our expression reduces to the familiar equilibrium result. In Sec. m we introduce the two 
classes of initial conditions that we will apply to unstable plasmas, and in Sec.jVjwe study the effect of self-interaction 
which needs to be subtracted from the energy loss formula. We develop our formalism in Sec. VI and in Sec. VII 
we apply it to extremely prolate and oblate systems. Our results are summarized and the outlook is discussed in 
Sec. |VIII| In Appendix we show that the energy loss is real, and the temporal axial gauge is compared with the 
Feynman-Lorentz gauge in Appendix [B| 


A preliminary account of our findings was published in the series of conference proceedings |2IH24| . In these reports, 
we used a specific choice of initial conditions and we were not aware that our results crucially depend on this choice, 
to the extent that the test parton can not only loose energy in an unstable plasma, but could also gain the energy, 
depending on how the initial conditions are chosen. Therefore, the results presented in do not reveal the full 

complexity of the problem. 

Throughout the paper we use natural units where h = c = 1. The indices i,j, k = 1, 2,3 and /i, = 0,1,2, 3 label, 

respectively, the Cartesian spatial coordinates and those of Minkowski space. 


II. GENERAL FORMULA 


Our formalism is based on the HL QCD effective action. It can be shown that the Wong equations uni and the 
linearized Yang-Mills (Maxwell) equations can be obtained directly from this action [TT] . The Wong equations describe 
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the motion of a classical parton moving through the fields of a plasma. The motion of the parton changes the field 
configurations, which is self-consistently taken into account through the linearized Yang-Mills equations which relate 
the chromodynamic fields to the parton charge and current. We emphasize that even though the Yang-Mills equations 
are linearized by the HL approximation, HL QCD is not equivalent to HL QED (up to an overall factor), because the 
gluons contribute to the color charge density and current in these equations. 

The Wong equations, which describe the motion of a parton in a chromodynamic field, are usually written in the 
Lorentz covariant form [10] 

u^(r), (1) 

(2) 

Qc(r), (3) 

where t, u^(t) and p^{t) are, respectively, the parton’s proper time, trajectory, four-velocity and four- 

momentum; and denote, respectively, the chromodynamic field strength tensor and four-potential in the 

adjoint representation of the SU(A^c) gauge group with the color index a = 1, 2,... N'^ — 1; g is the coupling constant, 
which is assumed to be small, and finally is the classical color charge of the parton. 

The Wong equations 0, § , and ([^ are supplemented by the linearized Yang-Mills equations describing the self- 
consistent generation of the chromodynamic field. We write the linearized Yang-Mills equations in a non-covariant 
three-vector notation where they have the familiar form of Maxwell equations in a medium. In Heaviside-Lorentz 
electromagnetic units, which are usually used in quantum field theory, we have 

V-Da(t, r) = pa{t,v) , V-Ba(t,r)=0, (4) 

VxE„((.r) = VxB.(*,r)=j.(t,r) + ?5^, (5) 

where E^, D^, are the chromoelectric field, chromoelectric induction and chromomagnetic field; Pa and are 
the density and current of the test parton, respectively. To close the system of Maxwell equations Q and ([^, the 
chromoelectric induction is expressed through the chromoelectric field 

Da{t,r) = J dt',v-r')El{t' ,y'), (6) 

where (t, r) is the chromodielectric permeability. We note that when the medium is on average color neutral, the 
chromodielectric permeability is proportional to <5“^, which can be factored out to produce an expression that carries 
no color indices as in Eq. 0. We also note that a quantity which depends on color indicies through 5°“^, or carries no 
color indices, is gauge independent. 

To solve the Wong equations 0, 0 , and 0 we adopt two simplifying assumptions. The first is that we choose 
the gauge condition 

u^,{T)Al^{x{T)) =Q, (7) 

which requires that the potential vanishes along the parton’s trajectory. Using this gauge, the third Wong equation 
0 simply states that the parton’s charge is a constant of motion, or that Qa is independent of r. The second 
important simplification comes from the fact that we consider a highly energetic parton and assume that its velocity 
V is constant and = 1. In an equilibrium plasma the characteristic momentum transfer |Ap| is of order gT and the 
parton’s momentum |p| ^ T, where T is the temperature. The hard loop approach requires gT T, and therefore 
|Av| X. |Ap|/|p| ^ 1. When we consider anisotropic systems, we assume the same hierarchy of scales which gives 
|Av| ^ 1. The physical picture is that due to interaction with the chromodynamic field the parton’s energy and 
momentum evolve in time without changing the magnitude of its velocity. 

Replacing the proper time r by the time t = yr and writing x''{t) = F{t) and = yf*, the first Wong equation 
0 gives r(t) = vt. Using this result, we obtain from the second Wong equation 0 with p = 0 

= 5Q“Ea(t,r(t)) • V. (8) 

Since the current generated by the moving parton equals 

ja(t,r) vt), 


dx^(r) 

dr 

dp^{T) 

dr 

dQajr) 

dr 


(9) 








we rewrite Eq. (|^ as 
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dE{t) 

dt 


/ 


d^rEa{t,r) • ja(t,r). 


( 10 ) 


To obtain the energy loss we must solve equations Q and ^ for the electric field and substitute into equation (101. 

The electric field that appears in equations (|^ and ( [To| ) is the total electric field, which is the sum of the external 
field generated directly by the moving test parton and the induced electric field produced by the charge distributions 
and currents that are induced by the parton in the plasma medium. The external electric field gives the parton’s 
self-interaction and does not contribute to the energy loss. The energy loss comes physically from the motion of 
the parton into the opposing induced electric field. We derive below an expression for the total electric field from 
Maxwell’s equations. At the end of the procedure, we must either show that the self-interaction does not contribute 
to the energy loss, or we must subtract it. 

It seems clear from Eq. (10) that if the parton moves into an electric field of opposite orientation to its current, the 
change in the energy will be negative and we have energy loss. We will show, however, that is not always the case. 
If the calculation is done as an initial value problem, then the sign of the energy transfer crucially depends on the 
choice of initial conditions. 


To solve Maxwell’s equations we use the usual method which is to Fourier transform the differential equations to 
change them into algebraic equations which can be easily solved. However, we do not use a standard (two-sided) 
Fourier transform. Our problem is to track the evolution of a parton starting from some arbitrary initial time (which 
we take to be t = 0) and calculate its behavior at future times. The non-equilibrated plasma is not time-translation 
invariant, and the energy loss formula should depend on the initial conditions, which means that we need to formulate 
the calculation as an initial value problem. In order to do this, we use a one-sided Fourier transformation defined as 


/(w,k) 



d3r-e*("’‘-k'-)/(t,r). 


/'00-\-icr 


' —oo+i<7 


duj 


d^k 

(27r)3 


e-*("’‘-k''-)/(a;,k). 


( 11 ) 

( 12 ) 


The inverse transformation (121 involves the real parameter cr > 0 which is chosen so that the integral over uj is taken 
along a straight line in the complex w-plane, parallel to the real axis and above all singularities of /(a;,k). 

The one-sided Fourier transform of the current 0 is obtained from Eq. 0 where the time integral is defined 
through the limit 


poo 

lim / 
£-*■0+ Jo 


i 

w — k • V -I- *0+ ’ 


(13) 


which gives 


ja(w,k) 


igQ°‘v 

w — k • V -I- f0+ ' 


(14) 


This procedure is mathematically equivalent to multiplying the current in equation (|^ by a factor which can 

be interpreted physically as imposing the boundary condition that the current goes to zero as the time approaches 
infinity. 

The one-sided Fourier transform of the relation 0 provides 

Dl{u;,k)=e^^{u:,k)Ei{u:,k). (15) 


Applying the one-sided Fourier transform 


to the Maxwell equations (|^ and ([^ and using the relation (151 gives 


iEe'‘^{oj,k)Ei{uj,k) = pa{uj,k), 


(16) 


iEB:{oj,k) = 0, (17) 

*e*^'=fc^A,^(a;,k) = (a;, k) + i?*Jk), (18) 

ze*^'=fc^H,^(co,k) = j:(a;,k) -zcoe*^ (cc,k)Fi(co,k) -iJ^Jk), (19) 
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where we have written i?Q^(k) = Bl{t = 0,k) and similarly for Z3g^(k). These initial values come from the time 
integrals of the time derivatives of fields after performing an integration by parts. The algebraic equations (16) - (19) 
are solved for the field i?*(w,k) 




( 20 ) 


where we have defined the matrix 


E*^(a;,k) = -k^5^^ + + a;V^(a;, k). 


( 21 ) 


In a quantum field theory formulation, one uses the retarded gluon polarization tensor n*'^(a;,k) instead of the 
dielectric tensor e^^iuj^k) and the two quantities are related to each other as 


e^^{uj,k) = - ^W^{uj,k). 




( 22 ) 


The full polarization tensor carries Lorentz indices (^, ^ = 0,1, 2, 3), which label coordinates in Minkowski space, and 
not Cartesian indices {i,j = 1,2,3). The components of the polarization tensor not determined by Eq. (22) can be 
reconstructed from the transversality condition fc^n^‘^(fc) = 0 with = (w, k) which is required by gauge invariance. 
Using the relation (22), the matrix (21) can be written 


(w, k) = <5*^ (w^ -k^) + Ek^ - (w, k) = (A-^ (w, k)) 


(23) 


This result is recognized as the inverse retarded gluon propagator in the temporal axial gauge (^° = 0), and we will 
write the propagator as A®l(a;,k). Although the matrix E has been derived with no reference to a gauge potential, 
the form of the gluon propagator depends on the chosen gauge. We return to this alleged conflict in Appendix 
where we also show that the temporal axial gauge is particularly convenient for our energy loss calculation because it 
naturally provides gauge independent results. To reach this goal in the Feynman-Lorentz gauge, which, in particular, 
was used in [7], current conservation must be explicitly enforced. 

The effects of the medium are contained in the dielectric tensor, or the polarization tensor. Performing a linear 
response analysis of kinetic equations in the collisionless limit, or equivalently working in the diagrammatic hard loop 
approximation, the dielectric tensor of a locally colorless anisotropic plasma equals [25l [26] 


e^^{uj,k) = 5^^ + 


(fp 


JL 

2uj J (27r)3 w — Vp • k + i0+ 


1 _ 

UJ ’ 


U) 


V«/(P), 


(24) 


where p and Vp = p/|p| are the momentum and velocity of a massless parton, and /(p) is the distribution function 
for hard partons in the plasma. For the SU(Ac) gauge group /(p) = n(p) + h(p) + 2Ncng{p), where n{p), h(p), 
ng(p) are the distribution functions of quarks, antiquarks and gluons of a single color component. We remind the 
reader that the chromodielectric tensor does not carry any color indices, as the state corresponding to the momentum 
distribution /(p) is assumed to be colorless. The prescription makes the Fourier transformed dielectric tensor 
e®-^(t, r) vanish for t < 0. In kinetic theory, the infinitesimal quantity i0“'" can be treated as a remnant of inter-particle 
collisions. Integrating by parts, the chromodielectric tensor (24) can be rewritten in the form 




f (Pp /(p) 


2a;^ J (27r)^ 






■ Vpk^ 


W - V„ 


k + i0+^ 


(k2 - 


w - v„ 


k + f0+)2j’ 


(25) 


which is often more convenient to use than the expression (24). 


The energy loss in equation (101 can now be written in terms of the Fourier transformed field and current. Performing 
the inverse transformation (12), we have 


dE{t) 

dt 




poo+io-' . 


d^k 


and substituting the formulas (14) and (20) into Eq. ( |26| ) , one obtains 

/■“+*'" doj dP 


dElt) 

dt 


' —OD+i<7 J —oo+ia' 


d^k 

(^ 




(26) 


(27) 






(k) 
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The integral over uj' can be done easily since the integrand has only one pole at uj' = —oj = —k • v. The result of 
integration over uj' is 


dE{t) 

dt 


gQ^v^ 


d^k 

(27r)3 


^oo+icr 
> — 00 -\-i (7 


27ri 


{E-y^{uJ,k) 


'iujgQ°‘v^ 

- UJ — UJ 




(28) 


which is the main result of this section. Equation (28) gives the change of energy of the parton as a function of time, 
and the expression depends on the initial conditions. The integral over uj is controlled by the poles of the matrix 
(or equivalently the gluon propagator) which determine the gluon collective modes in the system. These 
modes are found as solutions of the dispersion equation 


det[E(a;, k)] = 0. 


(29) 


Physically this means that the test parton does not interact with plasma constituents but rather with the plasma 
collective modes. 


In Sec. IV we discuss how to choose the initial conditions which enter the energy loss formula (28). In the next 
section we show that in the equilibrium limit Eq. (28) reduces to the familiar result which is independent of the initial 
conditions. 


III. EQUILIBRIUM LIMIT 


When the plasma is in equilibrium all collective modes are damped and all poles of the propagator Ab(cu,k) = 
jt) are located in the lower half-plane of complex uj. The corresponding contributions to the energy loss 


(28) exponentially decay in time, and the only stationary contribution is given by the pole w = w = k - v which comes 

which include the initial values of the 


from the current of the test parton. This means that the terms in Eq. 
fields can be neglected. It is mathematically equivalent to use a two-sided Eourier transform from the beginning of 
the calculation, which means that the initial conditions do not appear in the Maxwell equations ( |16[ ), and the Fourier 
transform of the current is just proportional to 5{uj — k • v). The result is that, once again, the only contribution 
to the integral over uj comes from w = w = k • v. In both approaches the result is that the energy loss of a high-energy 
parton traversing an equilibrium plasma is given by the time independent expression 


dE 

dt 




d^k 

(27r)3 


cc(E-i)*^-(a;,k). 


(30) 


Since the parton’s color charge is not an observable quantity because of its gauge dependence, the energy loss (30) 
has to be averaged over the parton’s color state. This is achieved by means of the relations 


dQQa = 0 , 


(31) 


and 


dQ QaQa — C 2 , 


(32) 


which are derived in m-, (72 = 1/2 for a quark in the fundamental representation of the SU(Vc) gauge group and 
C 2 = Nc for a gluon in the adjoint representation. Using the relation (32), the color averaged energy loss is 


dE ■ 2ry i i f 
- = -^gCnvv^J ^ 


cc(E-i)*^-(w,k), 


(33) 


where the color factor Cn, is given as 


Cr = 


^ for quark. 


Na 


C2=Nr. 


for gluon. 


It is easy to see that the result in equation (33) is real. Since the electric field and electric induction are both real 
in coordinate space, it follows from Eq. § that the dielectric tensor obeys the relations 


3fi£*^(-a;,-k) =3fi£*^(w,k). 


3£*'^(—w, —k) = —3 £*'^(w, k). 


(34) 
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Since the analogous relations hold for the matrix (propagator) S“^(a;,k), the real and imaginary contributions to 
the integrand in Eq. (331 are, respectively, odd and even as functions of k. Therefore, only the imaginary part of 
S“^(w,k), which is responsible for dissipative phenomena, contributes to the integral (33), and the energy loss is real 
as it should be. 

In an isotropic plasma the dielectric tensor can be decomposed into longitudinal and transverse components 


hi hj / 

£*'^(w,k) = eL{uj,'k) +er(w,k) 
and the matrix (w, k) can be inverted to obtain the propagator as 


- 


1^ 


(S-i)®^(a.,k) = 


1 1 
CLi2e^(a;,k) k^ a;2£T(w, k) — k^ 




Substituting this expression into Eq. (33), the energy loss is written 


H 


f d^k ui 

1 k^v^ - 

1 (27r)3 k2 

_£L(a;, k) a;2£2.(a;,k) - k2_ 


Using the symmetry relations (34) for £i,T(w,k), Eq. (33) becomes 

cPk I 




(27r)3 k^ 


Im£L(a;,k) w^(k^v^ — a;^)Im£T(w, k) 


|eL(w,k)|2 


+ 


|a;2£T(a;, k) — k^p 


(35) 


(36) 


(37) 


(38) 


As discussed under equation (10), the energy loss formula ( |10[ ), and consequently the formula (28), includes the self¬ 
interaction of the test parton with the electric field generated by the parton’s current The parton’s self-interaction 
should not contribute to the energy loss (37), and therefore we need to calculate this contribution separately and, if 
it is not zero, we need to subtract it from the energy loss obtained from Eq. (37). Since the effect of self-interaction is 
the same in a vacuum and in a medium, we derive it substituting into Eq. (37) the dielectric functions of the vacuum, 
which are 


£L(a;,k) = £T(w,k) = 1. 


(39) 


Using Eq. (39) the formula (37) gives 


dE 

dt 


- v^) 

.9^Cr 




d?k oj 

(27r)3 _ ]j2 

d{cos9) cos 9 


(27r) 


poo p-\-l 

(1 — v'^) / dkk 

Jo J-l 


I — cos2 


= 0 , 


(40) 


where we have chosen the axis ^ along the vector v and written a) = k • v = fc cos 9. Although the momentum integral 
is quadratically divergent, the angular integral vanishes and the three-dimensional integral is zero. The zero result 
is also expected from Eq. (38) because the vacuum dielectric functions (39) are purely real. Thus we see that the 
parton’s self-interaction does not contribute to the equilibrium energy loss formula (37) or (38). In Sec. |V]we will 


show that this is not the case when the energy loss calculation is formulated as an initial value problem. 


The result (37) or (38) agrees with the expression obtained in m using kinetic theory, and with the result for the 
energy loss due to soft collisions calculated in the HTL approximation [14] , see also the textbook m- However, this is 
not the complete energy loss but rather the soft contribution to it when the wave vector k is of the order of the Debye 
mass. Physically it corresponds to an interaction of the test parton with soft collective excitations of the plasma 
medium. The incompleteness of the formula (37) or (38) is signaled by the logarithmic divergence as |k| —>■ oo. To 
obtain the complete collisional energy loss, the formula (37) should be combined with the hard contribution describing 


elastic collisions of the test parton with plasma constituents with momentum transfer much exceeding the Debye mass. 
The hard contribution is not ultraviolet divergent, as the maximal momentum transfer is constrained by the collision 
kinematics. The soft contribution to the energy loss depends logarithmically on the upper cut-off k^^^ divided by the 
Debye mass m, while the hard contribution has a logarithmic dependence on the energy of the parton E divided by 
the same cut-off fcmax- The energy loss thus equals 


^^=Xln 

dt 


m 


Fin 


E 


(41) 




















































FIG. 1: (Color online) The parton energy loss per unit time in equilibrium plasma as a function of fcmax- 


It can be shown [HIlIl that the coefficients X, Y are equal to each other and therefore 


at \m 


(42) 


The result is that the cut-offs cancel and one obtains a good approximation to the energy loss from the soft contribution 
with the parton energy used as an upper cut-off. 

As was mentioned in the introduction, the energy loss in anisot rop ic QGP was computed previously by Romatschke 
and Strickland [S] . Their result can be obtained from our formula (33) by using an anisotropic propagator for (w, k) 
and including only the contribution from the pole to = uj. Clearly this procedure produces a result for the energy loss 
that is completely time independent. As we will see in the subsequent sections, the energy loss in anisotropic plasma 
is actually strongly time dependent because of the unstable modes. 

In order to compare our results for the energy loss in an unstable plasma to the corresponding equilibrium result, 
we have computed numerically the integral (37) in spherical coordinates. As already mentioned, the integral is 
logarithmically divergent at large k= |k|, so we introduce a cutoff k < When studying plasmas with massless 

constituents, the Debye mass m can be chosen as the only dimensionful parameter, and we therefore use a system of 
units where all dimensionful quantities are rescaled by the appropriate powers of m. Numerically one can simply set 
TO = 1. We define the Debye mass as 


2 _ 2 
m = g 


<Pp /(P) 

(27r)3 IpI 


(43) 


For isotropic momentum distributions, the formula ( |43[ ) coincides with the standard definition of the Debye mass. 
The same definition can also be used for the anisotropic plasmas studied in this paper. In Fig. we show the energy 
loss in isotropic QGP divided by g'^w? as a function of computed for Cr = Nc = 5 which corresponds to a gluon. 
Since the energy loss is divided by g'^w? we do not need to specify the value of g. The numbers from this figure will 
serve as a reference for our results on the energy loss in unstable plasmas. 


IV. INITIAL CONDITIONS 


When the plasma is anisotropic, the propagator A(a;,k) = E“^(a;,k) has poles in the upper half-plane of complex 
Lo which correspond to instabilities, and the contributions to the energy loss from these poles grow exponentially in 
time. This means that the terms in Eq. (28) which contain the initial values of fields Dq and Bq are amplified by an 


exponential factor and, in contrast to the equilibrium situation, they cannot, in general, be neglected. 
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A. Uncorrelated initial conditions 


The simplest choice of the initial condition is Dq = Bq = 0, which means that the energy loss formula (101 becomes 


dE{t) 

dt 




d^k 


f'OC-\-i(7 


' —oo+2cr 


doj 

27ri 


-e 




UJ 




(jJ — (jJ 


(44) 


where we have used the relation (321 to average over colors. In fact, the formula (44) holds for a whole class of initial 


conditions whenever Dq and Bg are independent of the test parton’s current. In this case, the contributions to the 
energy loss (28) which contain Dq and Bg are linear in the parton’s color charge Qa, and consequently they vanish 


when color averaging is performed using the relation (31). 


Physically this result can be understood as follows. Let us consider an electron moving in an external electromag¬ 
netic field which is independent of the current generated by the electron. The energy loss formula is given by the 
electromagnetic analog of the formula ^ where E(t,r(t)) is the external electric field along the electron’s trajectory. 
The electromagnetic analog of averaging over the parton’s color is the averaging over the possible charges of a hypo¬ 
thetical electron which could carry either negative or positive charge, or the averaging over the charges of an electron 
and a positron. If an electron’s energy increases by AE in the time interval At, a positron’s energy would decreases 
by —AE in the same field configuration and time interval. Therefore, after averaging over charges, the net change in 
the energy is zero. 


It is important to remember that the contribution to the energy loss from the first term in (28), which is proportional 


to the current and not the initial fields, is non zero, even when uncorrelated initial conditions are used. Mathematically, 
this happens because this term is proportional to the square of the charge. In an electromagnetic plasma, is strictly 
positive, and in a QCD plasma the factor does not give zero when averaged (see equation (32)). Physically we 


see that the energy losses of the electron and positron have same sign because they are not interacting with external 
fields which are independent of their currents, but instead with the electric fields which they have induced in the 
medium. We also note that the procedure of averaging over electric charges looks similar to that of averaging over 
colors but the physical situation is quite different. A color charge is gauge dependent and consequently it is not a 
physical observable. Therefore, the averaging over colors must be performed in order for the energy loss to have a 
physical meaning. 

It is interesting to note that we can obtain the same energy loss formula (44) in a different way. If we multiply 


the current in equation (|^ by a step function 0(t) and then repeat the whole calculation using the usual two-sided 
Fourier transformation, the identical result is found. The initial fields Dg and Bg do not appear in the two-sided 
Fourier transformed Maxwell equations, and the two-sided Fourier transform of the current with the additional step 
function is the same as the one-sided Fourier transform of original current. Although the same result can be obtained 
in two different ways, the physical interpretation of the two procedures is somewhat different. Using the two-sided 
transformation with the current multiplied by a step function, we assume that the plasma system exists for all times 
but the test parton appears in the plasma at t = 0. This was the problem studied in the papers liE]. When 
the one-sided Fourier transformation is used, it is understood that we observe the whole system, which includes the 
plasma and the test parton, starting only at f = 0. The initial values of the fields Dg and Bg can be chosen to be 
independent of the parton’s color state, but they could also be specihed differently. In the next section we consider a 
class of nontrivial initial conditions for which the helds Dg and Bg are strongly correlated with the current generated 
by the test parton. 


B. Correlated initial conditions 


We have shown in the previous section that if the initial conditions are chosen in any way that is independent of 
the parton’s current, they will not contribute to the energy loss. In this section we will consider another kind of 
initial conditions. First we note that although initial conditions are always required to solve differential equations, 
they are usually determined by physical arguments which go beyond the differential equations under consideration. 
We argue below that in the energy loss calculation, a kind of correlated initial conditions might be the most physical. 
We imagine that a process which is responsible for the occurrence of a test parton in a plasma system at the time 
t = 0, also polarizes the medium producing a chromodynamic held which is then correlated with the parton’s color 
state. We would like to see if the energy loss is sensitive to this kind of correlated initial condition. 

It is important to realize that the history of the system under consideration does not start at t = 0 when the test 
parton enters the plasma but earlier when the colliding nuclei begin to approach each other. When they collide, they 
produce among other particles the test parton. The plasma is produced at more or less the same time as the test 
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parton, and therefore, there is no reason to assume that the test parton and plasma, which are both a part of a 
bigger system, are completely uncorrelated. Note that this correlation is generated before the collision at t = 0, and 
therefore there is no violation of causality. 

Although the physical origin of a correlation between the parton and the plasma fields is clear, the question of how 
to express this correlation mathematically is much more difficult. To derive an upper limit, we assume that the parton 
enters the system in the remote past at t = —oo, observing that the parton’s current can be extended to the time 
interval from —oo to oo. Flying across the plasma, the parton polarizes the medium and induces a chromodynamic 
field. The initial fields Dq and Bg are identified with the induced fields at t = 0. 

To determine the fields Dg and Bg we solve the Maxwell equations and © using a normal (two-sided) Fourier 
transform with the time integral from —oo to oo. We use tildes to indicate that a two-sided Fourier transform was 
taken, which means that, for example, I?(a;,k) and Tl(a;, k) are different functions of uj but the same function of k. 
However, we note that e{uj, k) = e(w, k) and k) = k) because these functions obey the retarded initial 

condition and therefore e(t,r) = E“^(t,r) = 0 for t < 0. Solving the equations Q and © using a two-sided Fourier 
transform produces the result in equation (20) with and jl tilded and Hg^ = = 0, since the initial helds in 


Eq. (20) come from the t = 0 lower limit in the one-sided Fourier transform. Using the tilded version of the material 


relation (15), the electric induction is 


where the two-sided Fourier transform of the current in equation (|^ is 

ja(a;, k) = gQ°-v2TrS(uj - d)). 


Taking the inverse two-sided Fourier transform of the result (|45|), we obtain 

Dl 


/ OO 1 


(45) 


(46) 


(47) 


and setting t = 0, we arrive at 


DIM = -z5Q“d}e*^(d},k)(E-i)^'=(cc,k)u^ 
Using the same method we obtain the initial value of the chromomagnetic field 

BM) = -*5Q“e*'"fc'(S-')'='(cd,k)u'. 


(48) 


(49) 


The formulas (48) and (49) provide maximally correlated initial conditions. In order to consider initial conditions 
with differing degrees of correlation, we will multiply the initial helds (48) and (49) by a phase factor cos a S [—1,1]. 
The choices cos a = ±1 correspond to maximally correlated and anticorrelated initial helds. These two extreme cases 
provide limits on the possible effects of correlated initial conditions. We substitute the initial helds Dg and Bg given 
by Eqs. (48) and (49) into the energy loss formula (28) and insert the phase factor cos a as described above. After 
averaging over the parton’s color we obtain 


dEit) 

dt 


= ig^CRv'-k 


d?k 

(^ J- 






<50+2(7 


27rj 


+-^)*^(o.,k) 


(50) 




— cos a 


cu — cu 




This result, which reduces to the formula (44) when cos a = 0, will be further studied in the subsequent sections for 
two different unstable plasma systems. 


V. SELF-INTERACTIONS 


As already discussed in detail in the context of the equilibrium result (38), the energy loss formulas include the 
effect of self-interaction - also called the vacuum effect - which needs to be subtracted if it is non-zero. In this 
section we calculate th e sel f-interaction contribution to the energy loss given by Eq. (44) and ( [50| . We follow the 
same method as in Sec. Ill We evaluate the formulas (44) and (50) with the propagator S“^(a;,k) in the form (36) 
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with the vacuum dielectric functions (39). However, the calculation is not the same as the one done in Sec. Ill The 


equilibrium result (38) only has a contribution from the pole w = w, but the energy loss formulas (44) and (50) 


with a vacuum propagator also include contributions from the poles w = 0 and w = ±|k|, which make the effect of 
self-interaction time dependent. We discuss only the vacuum contribution to the energy loss formula for correlated 
initial conditions (501, because the corresponding result for uncorrelated initial conditions can be obtained by setting 
cos a = 0. 


To compute the vacuum effect we substitute into the formula (50) the vacuum propagator (36) with = 

k) = 1. In this way one finds the longitudinal part 


dExit) 


dt 


= ig^Cn 


d^k 


duj 


vacuum 


(27r)3 fc 2 2TTi 


and the transverse one 

dExit) 
dt 


= 


d^k uP- 
(27r)3 \~lP 


duj 

27ri w2 _ 


—oo+zcr 


LO 

UJ — UJ 


■ COS a 


cos a 


U)U1 + fc2 

0)2 — fc2 


(51) 


(52) 


where k = |k|. Performing the integral over w, which includes contributions from the poles at a; = a; and w = 0 in 
case of the longitudinal part (51), and the poles at w = w and w = Ek in case of the transverse one (521, we obtain 


dExit) 


dt 

dEx{t) 


dt 


= —{I — cosa)g^Cxi 
= — (1 —cosa) 


d^k u! sin(a;t) 


.. (27r)3 

ig^Cxi f d^k 


fc2 


('-S)( 


vacuum 


(27r)= 


k — UJ 


pi{k+uj)t 


(53) 


(54) 


We note that the pole uj = uj does not actually contribute to the transverse part (54) because the integrand is odd as 
a function of k and therefore it gives zero when integrated. 

The integrals over k are calculated in spherical coordinates using an upper cut-off fcmax- Summing the longitudinal 
and transverse parts, the complete vacuum contribution to the energy loss formula (50) equals 


dE{t) 


dt 


vacuum 


(1 - cos a)g'^Cn 


2(Si(/LniaxO Sin(fc,^ax0) E (2/Caaax^ Si(2/Caaax^)) 


where Si(z) is the sine integral defined as 


Si(,)= 

Jo a; 


(55) 


(56) 


The first term in the expression (551 is the longitudinal part and the second term represents the transverse piece which 
linearly diverges with increasing fcmax- Both the longitudinal and transverse contributions go to zero when t —>■ 0. 
The vacuum contribution to the energy loss formula with uncorrelated initial conditions (44) is given by Eq. (55) with 


cos a = 0. From Eq. (55) it is clear that the vacuum contribution is not zero unless we choose maximally correlated 


initial conditions (for which cos a = 1), and therefore the self-interaction effect must be subtracted from the energy 
loss formula in all other cases. 


VI. ANISOTROPIC PLASMAS 


The energy loss in isotropic plasmas has been calculated from the general formula (28) in section 


III In this section 


we develop our formalism to apply it to a general class of anisotropic momentum distributions of plasma constituents 
which was introduced in [ST] and has been used in various studies of QGP, see e.g. [32lI37) . These anisotropic 
distributions are obtained from the isotropic one by deforming it - squeezing or stretching - in one direction. The 
dispersion relations of the collective modes, which are needed to compute the energy loss, have been studied in 
great detail in our recent study [28] for all possible degrees of deformation from the extremely prolate case, when the 
momentum distribution is infinitely elongated in one direction, to the extremely oblate distribution, which is infinitely 
squeezed in one direction. 

In our derivation of the energy loss formula, and our calculation of energy loss for isotropic systems, we have mostly 
used the terminology of classical electrodynamics of continuous media, with the dielectric tensor playing a key role. 
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From now on we will switch to the language of quantum field theory, and make use of the polarization tensor and 
gluon propagator which were already introduced in Eqs. (22) and (23). The two languages are equivalent, as QCD in 
the hard-loop approximation is essentially classical, but the terminology of quantum field theory is more commonly 
used when working with the distribution introduced in |31j . 

In order to simplify the notation, in the rest of this section we omit the arguments which denote dependence on 
the wave vector. For example, we write a(w, k) as a{uj), A®'^ (a;,k) as A®'^ (w), etc. 


A. Propagator 


To compute the energy loss using the formula (28) we have to invert the matrix E defined by Eq. (21) or (23), which 


is the inverse gluon propagator in the temporal axial gauge. In isotropic plasmas the matrix depends on only one 
vector k. It can be decomposed into transverse and longitudinal components and is therefore easily inverted giving 
Eq. (36). We will now consider momentum distributions of the plasma constituents that can be obtained from the 


isotropic one by deforming it in one direction along the unit anisotropy vector n. In this case the matrix S depends 
on two vectors k and n, and it is symmetric = E-^*. To invert such a matrix, we introduce, following |38j . the 
vector riT 


nT = {S^ - -j^) 


(57) 


and we define four tensors m 


- 




= 


yy 


= 


I 


k 2 ’ 

= yn{ + yntp . 


The matrix E defined in Eq. ( |23[ ) is written 

E*^'(a;) = (uj^ - k^ - a{uj))A^^ + - -f{uj)C^^ - 5{w)D^\ 

where a, j3, 7 , 5 are the coefficients of the decomposition of the polarization tensor 

= a{w)A^^ + l3{uj)B^^ + + S{uj)D^^ . 

Writing the propagator A = E“^ in the same basis, the equation AE = 1 gives 
A*^'(a;) = Aa{uj){A^^-C^^) 

+ Ag(w) ((w" - k 2 - a{uj) - 7(w))S*2 + (^^2 _ , 

where 

A((^(a;) = — k^ — a(w), 

AQ^(a;) = —/3(a;)) — k^ — a(tLi) — 7 ( 0 ;)) — k^ny(5^(w). 


(58) 


(59) 


(60) 

(61) 

(62) 


From Eq. (60) we see that the poles of the propagator, which correspond to gluon collective modes, or plasmons, are 

A^i(a;)=0, Aci(a;)=0. (63) 


given by the dispersion equations 


The equations (63) are obviously equivalent to the general dispersion equation (29). 


B. Integrand 


Substituting the propagator (60) into the energy loss formula (50) and contracting all indices, we obtain an ex¬ 


pression that we will use to do calculations for the extremely prolate and extremely oblate momentum distributions 
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discussed in Sec. |VII| We use a spherical coordinate system with the z-axis along the anisotropy vector n. The angles 
9 and (p are the zenithal and azimuthal angles of the vector k, and 0 is the angle between the velocity of the test 
parton v and the anisotropy vector n. In our coordinate system the vectors n, v and k are 


n = (0,0,1), 

V = (sin 0,0, cos 0), 
k = A:(sin 0 cos^, sin 0 sine)), cos 0 ). 


(64) 


The energy loss formula (50) is written as 

dE{t) 


dt 


= 9^Cr 


d^k 


oo+ia- 7 

— Integrand, 


(27r)3 2TTi 


where the integrand is divided into several different pieces by writing it as 

Integrand = Aj H- Gj H- coscr^ylic Gic H” 


(65) 


( 66 ) 


The two terms Aj and Gj are the contributions from the first term in the square bracket in Eq. (28), which comes from 
the parton current (14|. When the initial fields are set to zero, or when we have uncorrelated initial conditions, these 
are the only terms that survive. They are proportional to A)(^(a;) and Ag^(w). The terms Gic, [^^]ic, [GG]ic 
are the contributions from the second two terms in the square bracket in Eq. (28) and come from the initial fields. 
They are proportional to Ayi(a;), AG(a;), Aa{oj) Aa{ui) and Ac{oj) Ac{ui). In the future we will refer to the first two 
terms in the formula (65) as ‘current contributions,’ and the last four terms will be called ‘field contributions.’ After 
performing all contractions we obtain 

Co\Cj'^{x^ — 1) — + l] 


A, = 


G, = 


{l-x^){A-Q)A-^\u) ’ 

Y‘^j3'{uj) — — l){(I3[/c^(a)^ ~ 1) ~ a{u}) — 7 ( 0 ;)] + 2kYS{uj)} 

{1 — x"^) {oj — oj)App^ {u}) 


-^ic — 

1—I 

I 

3^3 

Gic = 


III 

fc^(a) — w)(a) + (h) 

CjujA~^{uj) 

[GG]ic = 

„ , , kYC j il — X 


kYCj (1 — [P'{uj)S{uj) + P'{uj)8{uj)\ + Y"^ jd'{oj)(3'{uj) + k'^uj'^ (l — x'^Y 5(u))5{uj) 

w(l-a;2)Agi(a;)Agi(w) 


where we have used the symbols 

x = cos0, (xi = Lojk^ u} = uj/k^ Y = cos0 — xA, 
and defined the function /d'(u!) = uj^ — /3 (cj). 


VII. EXTREMELY PROLATE AND OBLATE PLASMAS 


In the early stages of a heavy ion collision, when partons are initially released from the incoming nucleons, the 
momentum distribution is strongly elongated along the beam - it has a prolate shape with the average transverse 
momentum much smaller than the average longitudinal one. Due to free streaming, see e.g. ISO], the distribution 
evolves in the local rest frame to a form which is squeezed along the beam - it has oblate shape with the average 
transverse momentum being much larger than the average longitudinal one. Prolate and oblate distributions can be 
obtained from an isotropic one by stretching or squeezing in the direction of a unit vector n which is chosen parallel 
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FIG. 2: (Color online) Unstable mode for the extremely prolate plasma: —a;?.(k) as a function of k and cos 6 in the domain 
where a;i(k) < 0. The angle 9 is between the vectors k and n. 


to the beam direction. In this paper we consider an extremely prolate and an extremely oblate distribution which are 
dehned as 


/ex-prolate (P) = HPt)— giPh), 

Pt 

/ex-oblate(p) = S{pL)h{pT), 


(67) 

( 68 ) 


where = p • n and Pt = |p ~ (p ' n)n|. The functions h{pT) and g{pL) are determined from the normalization 
condition (43) which now has the form 


2 9 

™ = tru 


g^ 

J d,pTh{pT) = ^j dpLgiPh)- 


(69) 


The collective modes are determined by substituting the distributions (671 and (681 into the dielectric tensor (24) 


which is related to the polarization tensor by Eq. (22 1 . One calculates the components of the polarization tensor (59) 


and solves the dispersion equations (63). The spectrum of plasmons of the extremely prolate and extremely oblate 
plasmas was analyzed in detail in our study |28] . Below we give a summary of the results. Solutions of the equation 


= 0 are called A-modes, 


and solutions of = 0 are referred to as G-modes. 


In the extremely prolate system, we can solve the dispersion equations analytically and write the propagators in 
the simple form 


,‘, 2^4 


Ag, (w) = 


(w"-w^)(A"-w^)(A"-wl), 


(70) 

(71) 


where the dispersion relations are 


= 


, m 


= X 


m 


2A± = 1 + a;^ ± \/(l — + 2m?/k"^). 


(72) 

(73) 

(74) 


The modes iwa, Au )2 and ±a;+ are pure real and exist for all wave vectors. The iw- modes are either pure real or 
pure imaginary. When k is greater than a threshold value, which we call /cpG (see Eq. (76)), they are real, and when 
k is less than this threshold they are imaginary and can be written as ±i^a with 7 ^ S M. 

Eor the oblate distribution, solutions of the dispersion equations can only be found numerically. There are 6 real 
modes (3 pairs) that exist for all wave vectors which we call iwa (A-modes), ±a;_ and ±w+ (G-modes). There are 
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FIG. 3: (Color online) Unstable modes for the extremely oblate plasma: —to'^ (panel (a)) and —(panel (b)) as functions of 
k and cos 9 in the domain where the modes exit. The angle 9 is between the vectors k and n. 


two thresholds (defined in Eq. @) below which imaginary modes appear. For k < fcoA there is a pair of imaginary 
A-modes and for k < koQ there is a pair of imaginary G-modes ±J 7 _. Note that we use the same terminology 

for prolate and oblate modes without introducing subscripts to distinguish them, but this will not cause confusion 


The threshold wave vectors for prolate and oblate systems are [3S] 


because the prolate and oblate systems are considered separately in sections |VIIA|and[VII^ 


kpG = ^|tan6>|, 


Tfi 

koA= -^\ cot ( 


(75) 


_m /1 cos0|\/cos^ 0 + 4 + cos^ 0 — 2 

^oG = ~o^\ ^ 2~^ ' 

2 V Sin 0 


When 0 —>■ 7 r /2 we have fcpc —>■ oo and the unstable prolate G-mode exists for all fc’s. This behavior is shown in Fig. 
for the prolate unstable G-mode. In the oblate system the situation is reversed and the threshold wave vectors fcoA 
and koG approach infinity when 0 —>■ 0. The unstable A and G-modes for the extremely oblate system are shown in 
Fig.0 

As discussed in Sec. [nj we calculate the frequency integral of the energy loss formula (65) with a contour in the 
lower half plane that encloses all singularities of the integral. The significance of the imaginary modes can be seen 
immediately. Denoting the real part of an imaginary mode generically as 7 , it is clear that the residue of a pure 
imaginary mode contains a factor e'*'* which grows exponentially with time. However, the magnitude of the unstable 
mode is small (in mass units), and the region of phase space for which the unstable mode exists is quite limited. Since 
there is also an oscillatory factor under the integral (28), it is not clear whether the energy loss will increase 

exponentially as a function of time. 

The integral over the wave vector k is taken numerically in spherical coordinates with the z-axis along the anisotropy 
vector n. Since the integral is ultraviolet divergent, we regulate it by introducing an upper cut-off at some finite 
momentum A:niax- For both oblate and prolate plasmas, there is a potential divergence when an imaginary mode goes 
to zero, as the wave vector approaches its threshold value. However, these divergences cancel exactly (sometimes in 
combination with the residue from the pole at w = 0). There are divergences that depend on the azimuthal angle 
when w = 0 and uj = ±a;_ but they are odd and can be regulated using a principal part prescription. 

The current contribution to the energy loss, or the energy loss with the uncorrelated initial condition, is very 
oscillatory and hard to calculate, but we have checked that it is of the order of the equilibrium energy loss discussed 
in Sec. |HI|and its magnitude is much smaller than the field contribution. One example of this current contribution is 
given in Sec. |VH A| 

In the two subsequent sections we present our numerical results on the energy loss in the extremely prolate and 
extremely oblate plasmas. In all our numeric al ca lculations Cn = 3, which corresponds to a gluon, and our results 
are expressed in the units of m. As in section HI the energy loss is divided by and therefore the value of the 

coupling constant g is not specified. 
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FIG. 4: (Color online) The integrand of the energy loss 
with cos a = 1 m extremely prolate plasma as a function 
of k and cos9 for 0 = 7r/12 and t = 8/m. The integral 
over azimuthal angle has been performed. 



as a function of t for 0 = 0, for two different choices of the fcmax- 
The red (solid) curve is for fcmax = 5m, and the blue (dashed) 
curve is fcmax = 3m. 


A. Extremely prolate plasma 


In Fig. I^we present the integrand of the energy loss in prolate plasma as a function of k and cos9 for 0 = 7r/12. 
The integral over azimuthal angle (j) has been done, and the small spikes at the top of the figure are produced by 
numerical issues. The meaning of the angles 9, </>, 0 is explained by Eq. ( [64| . Comparing this plot to that shown in 
Fig.d one clearly sees the influence of the unstable mode - the integrand is large in the domain of k and cos 9 where 
the mode co- exists. 

Fig. i shows the current contribution to the energy loss (or the energy loss with cos a = 0) as a function of time 
for 0 = 0. Since dE/dt is negative, the parton loses energy. The two curves represent two values of fcmax = 3m and 
^max = 5m, and one sees that the magnitude of the energy loss increases with k^^y^. The result is time dependent, 
but it is approximately of the same magnitude as the equilibrium energy loss at a given fcmax- As seen in Fig. the 
equilibrium energy loss equals — 0 . 125 ^m^ for fcmax = 3m and —O.lSg^m^ for fcmax = 5m. 

In Fig. [^we show the field contribution to the energy loss (with cos a = I) as a function of time for fc^ax = 5m 
and four angles 0 between the parton velocity and the anisotropy vector n. The energy loss dE/dt is positive and it 
increases exponentially with time, showing the effect of the unstable modes. The parton thus gains the energy and 
the magnitude of dE/dt at later times is much bigger than in equilibrium plasmas (see Fig. [^. The sign of the field 
contribution to the energy loss is determined by the sign of the phase factor cos a, and therefore if we change the 
initial condition from cos a = 1 to cos a = — 1 we will get exponentially growing energy loss instead of exponentially 
growing energy gain. Since the field contribution to the energy loss is much bigger than the current contribution, the 
sign of dE/dt is actually controlled by the sign of cos a. Therefore, the energy loss crucially depends on the initial 
condition. 


One observes in Fig. [^a strong directional dependence of the energy loss. For a prolate system, the most important 
wave vectors are those for which k T n, where the threshold wave vector (751 goes to infinity, see also Fig. When 
k T n the unstable mode has an associated electric field that is parallel to the vector n. This point is explained in 
Appendix The energy transfer is most efficient when the electric field is parallel to the velocity of the test parton 
(v II E). Therefore, one expects the largest energy transfer when v |j n. This argument is shown schematically in 
Fig. ii and verified by the results presented in Fig. which demonstrates that the magnitude of the energy loss is 
maximal at 0 = 0, and rapidly decays when the angle 0 grows. 

In Fig. I^we show (in a logarithmic scale) the energy loss as a function of fcmax for 0 = 0 and two times: t = 5/m 
= S/m. 


and t = 

as in the equilibrium case. As discussed in Sec. |III| 
theory. 


The energy loss oscillates slightly, but the fcmax dependence can be roughly approximated as logfcmax, 

the divergence at large fcmax indicates a breakdown of the classical 
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FIG. 6: (Color online) The field contribution to the energy 
loss in prolate plasma as a function of time for four angles 
between the parton velocity and the anisotropy direction: 
= 0 (red, dot-dashed), 0 = 7r/36 (orange, dotted), 0 = 
7r/12 (green, dashed), and 0 = 7r/6 (blue, solid). The cut-off 
parameter is fcmax = 5m. 



5 10 15 yt 20 


m 

FIG. 7: (Color online) The field contribution to the energy loss 
in prolate plasma at 0 = 0 as a function of fcmax for t = 5/m 
(circles) and t = 8/m (squares). The red (solid) curves are a 
logarithmic fit. 


B. Extremely oblate plasma 


Calculations in oblate plasma are much more difficult than those in prolate plasma because the components of 


the polarization tensor defined by Eq. (591 have a more complicated structure. They contain square roots that are 
not defined along the section of the real axis where the arguments of the roots are negative. There are therefore 
contributions to the frequency integral from the discontinuities between the upper and lower sides of the cuts that are 
difficult to calculate. We have checked for several cases that they are small when compared to the pole contributions 
and we therefore neglect them. One consequence of this more complicated structure is that the spectrum of collective 
modes is richer - there two unstable modes instead of one as in the case of prolate plasma. It is impossible to solve 
the dispersion equations analytically, and one can only obtain the dispersion relations numerically. Finally, there is 
a technical complication related to the fact that the dominant contribution to the energy loss in the oblate plasma 
comes from the domain of wave vectors k which are almost parallel to the anisotropy vector n. When k || n we have 
riT = 0 and the decomposition (58) is ill defined. This occurs because, when k||n, the matrix E does not depend on 


two independent vectors k and n but only on one vector k or n. Consequently, the decomposition (58) should include 
only two terms with the matrices A and B. The propagator has the form 


A*^(w) = 


1 


— k2 — a{uj) 




1 


w 


U}^ — P{uj) 


B^J 


(77) 


(a) (b) 




FIG. 8: (Color online) The configurations when the energy loss is maximal in the prolate (a) and oblate (b) plasmas. The 
arrows labeled with n, k, E and v denote, respectively, the anisotropy vector, the wave vector of the most important unstable 
mode, the electric field associated with this mode, and the velocity of the test parton. 




























18 



FIG. 9: (Color online) The field contribution to the energy loss in oblate plasma as a function of time for 0 = 7r/2. The red 
(dashed) line corresponds to the effect of ^-modes, blue (dotted) line is the G-modes, and the black (solid) line represents the 
sum. 


and the components of the polarization tensor are easily calculated in this special case as (see [5H1 for details): 






4a;4 


= 1 - ■ 


(78) 


In our calculation the domain where k||n was treated analytically and combined with the results of the numerical 
computation, as described below. Because of these technical difficulties, we give numerical results for the extremely 
oblate plasma only for a rather small value of fcmax = 2m. The equilibrium energy loss for this value of which 

will be used as a reference point, equals —0.079g^m^, see Fig.[^ 

As in the case of prolate plasma, the current contribution is significantly smaller than the field contribution. The 
latter is shown in Fig. when the parton’s momentum is perpendicular to the anisotropy vector n (0 = 7r/2). The 
red (dashed) line represents the contribution due to the A-modes, blue (dotted) line represents the G-modes and black 
(solid) one gives the sum. The black (solid) is not exactly the sum of the red (dashed) and blue (dotted) because in 
the calculation with all modes the points at x = ±1, which are obtained analytically, are combined and integrated 
together with the numerical data which is calculated over the range —0.9996 < x < 0.9996. One observes in Fig. 
that the unstable A-mode is responsible for the largest effect. Since the held contribution to the energy loss is much 
bigger than the current contribution, the sign of the energy loss is determined by the sign of cos a which expresses 
the dependence on the initial conditions. dE/dt is negative for cos a < 0 and it is positive when cos a > 0. As seen 
in Fig. the energy loss in oblate plasma can be orders of magnitude bigger than in an equilibrium plasma with the 
same fcmax- 


For an extremely oblate system, the most important wave vectors are those for which k || n, since both of the 
thresholds Icqa and k^G go to inhnity in this limit, see Eq. (76). This behavior is also shown in Fig[^ As explained in 
Appendix instead of two different pairs of imaginary modes A and G, we have two pairs of identical modes which 
are purely transverse when k || n. The electric field associated with these modes is perpendicular to both k and n. 
Since the energy loss is maximal when the parton velocity is parallel to the electric field, such a situation occurs in 
the oblate system when v is perpendicular to n, or 0 = 7r/2. This argument is shown schematically in Fig. [^. The 
effect is seen explicitly in Figs. [T0| and 13 The left panel of Fig. [I0| shows that for both A- and G-modes the energy 
loss is dominated by the region a: « 1, and the right panel proves that when x = \ the biggest effect is observed when 


0 = 7r/2. Fig. 11 presents the energy loss as a function of 0 for t = 25/m. The figure shows that dE/dt drops rapidly 


when 0 becomes smaller than 7r/2. 


VIII. DISCUSSION AND OUTLOOK 


Let us first summarize our study. We have derived the energy loss formula for a high-energy parton flying across an 
unstable plasma which experiences a rapid temporal evolution due to exponentially growing collective modes. Since 
the formula is found as the solution of an initial value problem, initial values of the chromodynamic fields present in 
the plasma must be chosen. Except in special cases, the energy loss formula includes an effect of self-interaction which 
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FIG. 10: (Color online) The integrand of the field contribution to oblate energy loss after integrating over (j> and k for t = 25/m. 
In the panel (a) the integrand is shown as a function of a; = cos6' for 0 = 7r/2 and in the panel (b) as a function of 0 for a; = 1. 
The red (dashed), blue (dotted) and black (solid) lines in the panel (a) represent the effect of, respectively, A-modes, G-modes 
and the sum of A plus G-modes. 


must be subtracted to get a physically meaningful result. In the case of equilibrium plasmas, the initial conditions 
are ‘forgotten,’ and the well-known formula of collisional energy loss is reproduced. When the initial conditions are 
chosen in such a way that the initial fields are not correlated with the current generated by the test parton, the parton 
typically looses energy, and the magnitude of the energy loss is comparable to that in an equilibrium plasma of the 
same mass m (431. When the initial chromodynamic field is induced by the parton, it can be either accelerated or 
decelerated depending on the relative phase factor. With correlated initial conditions, the magnitude of the energy 
loss grows exponentially in time and can much exceed the absolute value of the energy loss in an equilibrium plasma. 


We have derived an expression for the energy loss for arbitrarily prolate or oblate plasmas, and performed numerical 
calculations for the specific examples of the extremely prolate and extremely oblate systems. The energy loss is not 
only time dependent but it is also strongly directionally dependent. The configurations when the energy loss is 
maximal in the prolate and oblate plasmas are illustrated in Fig. In these special configurations, the magnitude 
of the energy loss can be much bigger than that in an equilibrium plasma. Beyond a narrow cone which is centered 
around the optimal direction, the energy loss rapidly drops. 


It is interesting to consider the possible consequences of our findings for the jet suppression observed in relativistic 
heavy-ion collisions. Since a high-energy parton can be either accelerated or decelerated in an unstable plasma, we 
expect that the energy loss strongly fluctuates, and that the fluctuations are particularly large in the configurations 
depicted in Fig. Quark-gluon plasma at an early stage of a relativistic heavy-ion collision has initially a prolate 



FIG. 11: (Color online) The field contribution to the energy loss in oblate plasma as a function of the angle 0 for t = 25/m. 
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momentum distribution which evolves fast due to free streaming to an oblate momentum distribution. During the 
process of equilibration the plasma is oblate and it remains oblate in the subsequent evolution because of viscosity 
effects [10]. Jet quenching is observed at both RHIC and LHC at almost vanishing rapidity in the center of mass of 
colliding nuclei. This configuration is just as shown in Fig. j^, where the jet momentum is transverse to the vector n. 
We suspect that the jet quenching pattern can be changed when the jet axis is tilted in such a way that the near-side 
jet has a small but positive (negative) rapidity while the away-side jet has a small but negative (positive) rapidity. 
The effect of unstable modes is then reduced and the energy-loss fluctuations are expected be much smaller. 

One should remember that we have discussed here only collisional energy loss. There are simple arguments that 
indicate that radiative energy loss behaves very differently |20j . Radiative energy loss is controlled by the parameter q 
which measures the momentum broadening of a parton. This parameter is by definition positive and grows exponen¬ 
tially in an unstable plasma, as does the radiative energy loss, which is always negative. Therefore, before we draw a 
conclusion about the possible role of unstable plasma in jet suppression phenomenology, the effects of both collisional 
and radiative energy loss must be combined. This requires a computation of q in unstable plasmas which we plan to 
publish soon. 
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Appendix A: Reality of energy loss 


We prove here that the energy loss (28) is real for any momentum distribution that satisfies the mirror symmetry 


/(p) = /(~P)- For this purpose we take the complex conjugate of the formula (28) and obtain 


dE*{t) 

dt 


= / tS (a;, k) [ - + e^'='fc'=<(k) - a;D*^(k)], (Al) 


where the inverse matrix E ^ is replaced by the retarded propagator A. Now we change the integration variables 
Lo —^ —u) and k —^ —k which gives 


dE*{t) 

dt 


= -gQ^v^ 


d^k 

(27r)3 


nOC-\-i(T 


> —oo+2cr 


duj 

‘I'Kl 


-*(“-“)*A”^'(-a;, -k) [ - + ujDl^i-k)] . (A2) 


Since the initial fields Boa(i") and Doo(r) are pure real in coordinate space, we have 

Boa(k) = B*J-k), Doa(k) = D*J-k). 


(A3) 


In our study |28j we have proven that for mirror-symmetric momentum distributions, the retarded propagator defined 
by Eq. (23) satisfies the relations 


5RA*^(-w,-k) = 3fiA*J(a;,k), 


SA*J'(-a;,-k) = -S5A*^(w,k), 


which give 


A“^'(-a;,-k) = A*^(a;,k). 


(A4) 


the proof that the energy loss given by the formula (28) is real. 


Using the relations (A3) and (A4), the right side of Eq. (|A2|) is identical to the right side of Eq. (28), which completes 


Appendix B: Temporal axial and Feynman-Lorentz gauges 

In this appendix we show that the temporal axial gauge is particularly convenient for the energy loss calculation 
because it naturally provides a gauge independent expression for the energy loss that depends only on the electric and 
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magnetic fields. In contrast, current conservation must be explicitly enforced in Feynman-Lorentz gauge. To simplify 
the problem, we consider here an electromagnetic plasma. 

The electromagnetic analog of the energy loss formula ( [2^ is clearly gauge invariant, as the derivation of the 
formula is gauge invariant at every step. The gauge dependent potential is not used at all, and the energy loss 
is written in a form that depends only on the dielectric tensor and electric and magnetic fields, which are physical 
quantities. However, when we switch to the terminology of quantum field theory and the inverse dielectric tensor is 
replaced by the photon propagator in the temporal axial gauge, the gauge independence of the formula (28) is not 


evident any more. In this section we will show that although the energy loss formula looks different in the Feynman- 
Lorentz gauge, it is still gauge invariant. We also explain why temporal axial gauge is much more convenient for the 
energy loss calculation. In this appendix we use the usual (two-sided) Fourier transformation and not the one-sided 
transformation, which was used in Sec. [H] 

To further simplify the problem we will consider not the whole energy loss formula but only the electric field 
generated by the test particle in vacuum. We will solve the Maxwell equation 


d^F>^‘'ix) = nx), 


(Bl) 


where x = {t, r), — d'^A^ and is the particle’s current. The electric field, which is the physical quantity 

of interest, is expressed through the four-potential as 


E(a;) = -WA°{x) - A{x), 

which in momentum space is 

E{k) = -ikA°{k)+iujA{k). 

We note that k = (a;,k) denotes here the four-vector and not |k|. 


(B2) 


(B3) 


In order to solve equation (BI) for the potential, one must choose a gauge. The resulting solution has the form A 


= propagator x current. Both the propagator and the vector potential are gauge dependent. However, if we calculate 
the electric field from the potential using Eq. (B2| or (B3), the result must be gauge independent. This is true when 
current conservation is imposed. 


We start by considering Feynman-Lorentz gauge [d^A^ = 0) in which the Maxwell equation (Bl) is 

□B^(x)=r(x), 

and the (two-sided) Fourier transformed solution reads 

A^^ik) = A^/^^ik)Mk), 

where 




g‘ 




k^ + isgn(w)0+ 


= g^'^DFLoik) 


(B4) 


(B5) 


(B6) 


is the retarded photon propagator in the Feynman-Lorentz gauge. From equations (B3), (B5) and (B6) we obtain the 
electric field generated by the current j‘'{k) 


E\k) = -iDFFG{k){Fj\k)-wj\k)). 


(B7) 


Now we consider the temporal axial gauge (B° = 0). The (two-sided) Fourier transformed field equation (Bl) splits 
into two equations 


The solution of the second equation 


IS 


where 


^tag(^) ~ 


1 




uF F isgn(a;)0+ k^ 


-ujFA^k) = f{k), 

(B8) 

-Fk^jA^ik) = j\k). 

(B9) 

■^TAg(^)7^(^)’ 

(BIO) 

- -h ^ kk^\ 

+ isgn(w)0+ V k^ / 

(Bll) 
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is the retarded photon propagator in the temporal axial gauge, 
using = 0) we obtain 


Substituting the solution (BIO) into Eq. (B3) and 


E"{k) = (k). 


(B12) 


We have found that in Feynman-Lorentz gauge the electric field is given by Eq. ([B7j) and in temporal axial gauge 
it is given by Eq. (B12) with the additional constraint (B8). The two equations for the electric field look different, 


but if current conservation is imposed they are in fact the same. Current conservation gives the relation 


ujf{k) = k • j(fc). 


(B13) 


Using (B13) the electric field obtained from the Feynman-Lorentz gauge (B7| can be written in the form 


E\k) = -DFLG(fc)(fc*fc^' - u}‘^S^^)fik). 

UJ 


(B14) 


Equations (B6) and (BID give the equality 


A¥AG(fc) = ^^FLG(fc)(fc*fc^' - 


(B15) 


Using Eq. (B15l it is easy to see that the expressions (B12) and (B14| are equivalent. 


When working in the temporal axial gauge, current conservation merely tells us that the solution ( B10[ ) satisfies 
Eq. (B8) automatically, and the electric field is naturally gauge independent. Equivalently, the electric field in (B12) 
can be derived from the Maxwell equations (|^ and ([^ without any reference to the four-potential . In contrast, 
if Feynman-Lorentz gauge is used, current conservation must be explicitly enforced. The authors of [7] resolved this 
problem by modifying somewhat artificially the parton’s current. 


Appendix C: Important configurations 


In this appendix we look at the prolate system in the special case that the wave vector of the unstable mode is 
perpendicular to the direction of anisotropy (k T n), and the oblate system when these two vectors are parallel 
(k II n). These regions of k are important because they are the part of the domain of k for which the unstable modes 
exist up to infinite k, see Eqs. (75) and (76) and Figs. [^and[^ We will further show that the energy loss is maximal 
when the velocity of the test parton v is parallel to n in the prolate plasma and when v T n in the oblate one. The 
arguments discussed in this appendix are illustrated in Fig. 

We start with the prolate system. The linearized Yang-Mills or Maxwell equations of electric field E(w,k) can be 
written as [25] 


E*^'(a;,k)U^'(w,k) = 0, 


(Cl) 


with the matrix E defined by Eq. (21). Since the equation 
of the matrix E vanishes 


is homogeneous, there are solutions if the determinant 
this is the general dispersion equation (|^. When n = (0, 0,1) and k = {k, 0,0), Eq. (Cl) 


IS 


where 


— l3{uj) 0 

0 — k'^ — a{uj) 

0 0 


a{uj) = P{ui) 



0 
0 

Q;(a;) — 7 ( 0 ;) 


7(w) = 



1- 

1_ 


Ey(uj,k) 

- 

_ U^(a;, k) _ 




= 0 , 


2^2 


(C2) 


(C3) 


The imaginary modes appear as solutions of the equation — aiio) — j(uj) = 0 which controls the z-component 

of the electric field. Therefore, the exponentially growing component of E is parallel to n. Since the maximal energy 
loss occurs when v |j E, the maximal effect requires v |j n. 
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Let us now consider the oblate plasma. When n = (0,0,1) and k = (0, 0, k), Eq. (|C1|) is 


- 1 

e 

to 

1 

to 

1 

Q 

o 

o 


1 

(S' 

H 

0 — a{uj) 0 


Ey{u;,'k) 

1 

'll 

1 

CN 

3 

o 

o 

_1 


- 1 

_1 


(C4) 


where the coefficients a(w) and /3(w) are given by Eq. (78). The imaginary modes appear as solutions of the equation 
— C({uj)) = 0 which controls the x- and ^-components of the electric field. Therefore, the exponentially 
growing component of E is perpendicular to n and the maximal energy loss occurs when v T n. 
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